Metastable States in CA models for Traffic Flow 
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Abstract. Measurements of traffic flow show the existence of metastable 
states of very high throughput. These observations cannot be reproduced 
by the CA model of Nagel and Schreckenberg (NaSch model), not even 
qualitatively. Here we present two variants on the NaSch model with 
modified acceleration rules ('slow-to-start' rules). Although these mod- 
els are still discrete in time and space, different types of metastable states 
can be observed. 



1 Introduction 

Simulation results of the CA model for traffic flow introduced by Nagel and 
Schreckenberg [0 (NaSch model, for a detailed explanation see show a 

reasonable agreement with experimental data, although the model definition 
includes only a few simple update rules. Due to its simplicity it can be used very 
efficiently for computer simulations, but also an analytical treatment is possible 

One important experimental result pj which cannot be gained from the 
NaSch model is the occurance of metastable states near the density of maxi- 
mum flow. Very recently such states have been found in a modified version of 
the NaSch model where a limited braking capacity of the cars and a continuous 
space is considered pj . The continuous space allows very small amplitudes of the 
fluctuations compared to the discrete NaSch model. Therefore one could ask if 
a continuous space is neccessary to obtain metastable states. 

Here we present the results of a numerical investigation of two variants on 
the NaSch model with modified acceleration rules, which are still discrete in time 
and space. First, we assigned a larger braking probability to stopped cars and 
second, we analysed the modification suggested by Takayasu and Takayasu B. 
It turns out that metastable states can be found in both cases. 



2 NaSch model with 'slow-to-start' rule 



In this part of the investigation we consider the NaSch model with maximum 
velocity v max — 5 and braking probability p = 0.01 of the moving cars. In 



contrast to the original model, a higher value p s = 0.5 of the braking probability 
is assigned to stopped cars. 

We performed simulation runs for periodic systems containing L = 10000 
lattice sites. This system size is sufficient to exclude finite size effects for the 
chosen set of parameters. Fig. p] shows the fundamental diagram of the modified 
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Fig. 1. Fundamental diagram of the modified NaSch model (v max = 5, p s = 0.5, 
p = 0.01, L = 10000). 

model. Obviously the average flow J(p) can take two values in the density interval 
between p\ and pi dependent on the chosen initialisation. The larger values of the 
average flow can be obtained using a homogeneous initialisation of the system. 
The lower branch is obtained starting from a complete jammed state. Moreover 
varying the particle number adiabatically, one can trace a hysteresis loop. One 
gets the upper branch by adding cars to the stationary state with p < pi and 
the lower one by removing cars from the stationary state with p > p2- For a 
fixed value of p, AJ = J(p2) — J(pi) depends linearly on p s for wide a range of 
parameters. 

The microscopic structure of the jammed states in the modified model differs 
from those found in the NaSch model. While jammed states in the NaSch model 
contain with an exponential size-distribution one can find phase separation 
in the modified model. The reason for this behaviour is the reduction of the 
outflow from a jam. Therefore the density in the free flow regime is smaller 
than the density of maximum flow and cars can propagate freely in the low 
density part of the lattice. Due to the reduction of the density in the free flow 
regime, no spontaneous formation of jams is observable in the stationary state, 
if fluctuations in the free flow regime are rare. 

This picture can be confirmed by a simple phenomenological approach. Ob- 
viously the flow in the homogeneous branch is given by Jh = p(v max —p) = pvf, 
because every car can move with its desired velocity v / . Assuming that the high 
density states are phase separated, we can obtain the second branch of the fun- 



damental diagram. The phase separated states consist of a large jam and a free 
flow regime, where each car moves with velocity Vf. The density in the free flow 
regime pf is determined by the average waiting time T w = jj of the first 
car in the jam and Vf, because neglecting interactions between cars the average 
distance of two consecutive cars is given by Ax — T w Vf + f = pj 1 . Using the 
normalisation L = Nj + NpAx (Nptn number of cars in the free flow regime 
(jam)) we find that the flow is given by 

J s (p) = (l-p)(l-p s ). (1) 

Obviously pf is precisely the branching density pi, because for densities below 
pj the jamlcngth is zero. It should be noted that this approach is only valid if 
{jps 3> p) and v max > 1 holds. Small values of p have to be considered in order 
to avoid interaction between cars due to the velocity fluctuations and v ma x > 1, 
because otherwise cars can stop spontaneous in the free flow regime and therefore 
initiate a jam. 




Fig. 2. Time-dependence of the jamlength for p = 0.095. The left part of the 
figure shows the time evolution of the length L(t) of one sample. The average 
over 10000 samples (right part of the figure) shows that the jamlength decays 
exponentially. 

Measurements of the average flow show that the lower branch of the fun- 
damental diagram is not stable near the density p\. Therefore we performed a 
more detailed stability analysis of the homogeneous and the jammed state near 
pi and pi. Near p\ the large jam present in the initial configuration resolves, 
the average length decays exponentially in time (Fig. ||) . It should be noted that 
this behaviour is not the consequence of a continuous "melting" of the large 
jam. In contrast, the jamlength is strongly fluctuating without any systematic 
time-dependence. Once a homogeneous state without a jammed car is reached, 
no jam appears again. Therefore the homogeneous state is stable near p\. 
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Fig. 3. Space-time diagram of the modified NaSch model for p — 0.15, L — 
400, p = 0.01 and p s = 0.5. The homogeneous initial state is not destroyed im- 
mediately, but after approximately 90000 lattice updates. In the outflow regime 
of the jam the density is reduced compared to the average density. 

Analogous to the metastable jammed states near pi, homogeneous initialisa- 
tions for densities slightly above p2 lead to metastable homogeneous states with 
short lifetimes. Fig. |3| shows the spontaneous formation of jams due to velocity 
fluctuations. The finite lifetimes of the homogeneous states are the qualitative 
difference between this model and the cruise-control limit [^J, where the time 
evolution of homogeneous states at low densities is completely deterministic. 

3 T 2 model 

Takayasu and Takayasu (T 2 ) H suggested a cellular automaton model with 
another type of a slow-to-start rule, which is defined as follows: A standing car 
with exactly one empty cell in front of it accelerates with probability q t — l~p t , 
while all other cars accelerate deterministicly. The other update rules of the 
NaSch model are unchanged. Due to this modification already for v max — 1 the 
particle-hole symmetry is broken. In Fig. ^ we show the fundamental diagram 
for v m ax = 1, p = 0.5 and pt = 1, i.e. stopped cars can only move if there are 
at least two empty cells in front of it. Obviously completely blocked states exist 
for densities p > 0.5, where the number of empty cells in front of each car is 
smaller than two. Since fluctuations are absent in those states, they have an 
infinite lifetime. Therefore the flow in the stationary state is zero. 

Although the states with a finite flow are not stationary, typically one has 
to perform an extremely large number of update steps until the flow vanishes 
for large system sizes and densities slightly above p = 0.5, because the number 
of blocked configurations is very small compared to the total number of config- 
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Fig. 4. Fundamental diagram of the T 2 model (v max = 1, p = 0.5, p t — 1) for 
two system sizes. For densities slightly above p — 0.5 the stationary state could 
only be reached for the small system. 

urations. Precisely at p — 0.5, the blocked state is unique and the typical time 
to reach this state diverges exponentially with the system size. Therefore we 
used very small systems in order to obtain the lower branch of the fundamental 
diagram. 

Another interesting feature is the form of the fundamental diagram. There 
is some experimental evidence that in certain situations the shape of the fun- 
damental diagram differs from the standard convex form. This behaviour of the 
average flow can be easily obtained tuning the parameter pt (2||7| . 

4 Summary 

The numerical analysis of two modifications of the NaSch model shows that 
metastable states can also be found in CA models for traffic flow. 

The NaSch model with 'slow-to-start' rule shows the coexistence of phase 
separated and homogeneous states in a density interval near the density p-i of 
maximum flow. Near pi interactions between cars become important and one 
can find spontaneous formation of jams. Therefore the reduction of the density 
in the outflow regime of a jam leads to stable phase separated states. The re- 
duction of interactions between cars in the free flow regime can be confirmed by 
a phenomenological approach, which gives very accurate results. In contrast to 
the cruise-control limit of the NaSch model j|] , where also metastable states can 
be found, fluctuations are present in both coexisting states. Finally we mention 
that this modified version of the NaSch model can reproduce the fundamental 
diagram of the spacecontinuous approach by KrauB et. al. |4j. 

In the pt = 1 limit of the T 2 model we find metastable states for densities 
beyond p = 0.5. For p > 0.5 states with a finite flow have very large lifetimes, but 
if the system reaches a blocked state these states are stable because of the absence 



of fluctuations. Therefore the pt = 1 limit of the T 2 model is complementary 
to the cruise-control limit of the NaSch model, where one finds the absence of 
fluctuations for homogeneous states at low densities. 

These modifications of the NaSch model show that CA models for traffic can 
reproduce a broad spectrum of experimental results although the update rules 
oversimplify the individual behaviour of the drivers. 
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